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Abstract 



For electronic excitations in the ultraviolet and visible range of the electromag- 
netic spectrum, the intensities are usually calculated within the dipole approxima- 
tion, which assumes that the oscillating electric field is constant over the length 
scale of the transition. For the short wavelengths used in hard X-ray spectroscopy, 
the dipole approximation may not be adequate. In particular, for metal K-edge 
X-ray absorption spectroscopy (XAS), it becomes necessary to include higher-order 
contributions. In quantum-chemical approaches to X-ray spectroscopy, these so- 
called quadrupole intensities have so far been calculated by including contributions 
depending on the square of the electric-quadrupole and magnetic-dipole transition 
moments. However, the resulting quadrupole intensities depend on the choice of 
the origin of the coordinate system. Here, we show that for obtaining an origin- 
independent theory, one has to include all contributions that are of the same order in 
the wave vector consistently. This leads to two additional contributions depending 
on products of the electric-dipole andelectric-octupole and of the electric-dipole and 
magnetic-quadrupole transition moments, respectively. We have implemented such 
an origin-independent calculation of quadrupole intensities in XAS within time- 
dependent density-functional theory, and demonstrate its usefulness for the calcu- 
lation of metal and ligand K-edge XAS spectra of transition metal complexes. 
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1 Introduction 



X-ray spectroscopy [T143] is a powerful spectroscopic tool for the elucidation of structural 
and electronic properties of materials |4j^6] and (bio-)molecular systems [7[j9]. X-ray 
absorption spectroscopy (XAS) probes the excitation of core electrons. Here, one dis- 
tinguishes excitations to low- lying unoccupied states (so-called prepeaks), excitations to 
states close to the ionization threshold (X-ray absorption near-edge structure, XANES), 
and excitations to continuum states (extended x-ray absorption fine structure, EXAFS). 
On the other hand, X-ray emission spectroscopy (XES) considers the emission of X-ray 
radiation after the formation of a core hole. 

Of particular interest are the applications of these techniques to study catalytic mech- 
anisms in situ (for examples, see, Ref. [3j [7j [l0| [IT] ) and to investigate biological or 



biomimetic systems containing transition metal centers (see, e.g., Refs. 12 18]). Usually, 
XAS and XES are used as fingerprint techniques in such studies to identify the oxidation 
state, spin state, and local coordination environment of a metal center. This requires 
either the comparison with spectra measured for model compounds or with theoretical 
predictions. To be able to extract additional information from X-ray spectroscopic mea- 
surements, the development of theoretical methods for the calculation of X-ray spectra is 



essential 19 



For EXAFS spectra, approaches based on scattering theory are well established and make 
it possible to extract structural parameters such as distances and coordination num- 



bers 20 21 . In contrast, for describing prepeaks and the XANES region in XAS spec- 
troscopy and for predicting XES spectra, quantum-chemical approaches are usually re- 
quired. To this end, a wide range of quantum-chemical methods have been developed for 
describing excitations from core orbitals (for reviews, see, e.g., Refs [l9]|22j[23]). Widely 



used are the static-exchange approximation (STEX) 26 , approaches based on tran- 
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sition potential density-functional theory (DFT) [27 28 , and time-dependent density- 
functional theory (TD-DFT) (23]. Within TD-DFT, core excitations are accessible by 
restricting the space of occupied-virtual orbital excitations (restricted-channel approxima- 



tion) Y29 30 , by selectively targeting excitations within a specific energy window 31 34 



by using a complex polarization propagator 1 35 , 36 , or with real-time TD-DFT meth- 
ods |37) . Recently, coupled-cluster response theory has also been extended to X-ray 
spectroscopy 



38.39 



By combining X-ray spectroscopy with quantum-chemical calculations, it becomes possi- 
ble to extract information on the electronic structure of molecular systems. For instance, 
the analysis of the prepeak intensities in ligand K-edge XAS spectra of transition metal 
complexes (i.e., excitations from the ligand Is to metal d orbitals) provides insights into 



the covalent contributions to metal-ligand bonding 40 46 . Another example is metal 



K-edge XAS, probing excitations from metal Is to d orbitals, which can be used to assign 



coordination numbers 47-49 and to probe details of the metal-ligand bonding mecha- 
nisms [l8|[44|[50| . Recently, we have demonstrated that the prepeaks in Fe K-edge XAS 
spectra of ferrocene derivatives are sensitive to subtle differences in the electronic structure 
at the iron atom, which are induced by substituents at the cyclopentadienyl rings, i.e., 
beyond the first coordination shell of the metal center [51 1. Such studies are facilitated by 
high-energy resolution fluorescence detection (HERFD) techniques, which can resolve the 



prepeaks with a much higher resolution than conventional XAS measurements 52 -54 



However, beside the challenges posed by the accurate quantum-chemical prediction of the 
absolute or relative energies of core excitations, for the prepeaks in K-edge XAS spectra, 
theoretical X-ray spectroscopy faces an additional problem. For the calculation of XAS 
intensities for hard X-rays, the well-known dipole approximation, in which the oscillator 
strengths are proportional to the square of the electric-dipole transition moments is not 
sufficient. The dipole approximation is based on the assumption that the wavelength 
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of the electromagnetic radiation is large compared to the size of the core orbital. For 
the high-energy radiation used in hard X-ray spectroscopy, this is not the case anymore. 
This is particularly important for the prepeaks in metal K-edge XAS spectra of transition 
metal complexes, which are usually dipole forbidden or have a very low dipole intensity 
due to symmetry. Thus, the intensity of these prepeaks is due to contributions that are 



not included in the dipole approximation 55 57 



Currently, contributions to XAS or XES intensities beyond the dipole approximation, so- 
called quadrupole intensities, are calculated by including additional contributions that are 
proportional to the squares of the electric-quadrupole and the magnetic-dipole transition 



moments 58 59 . However, these additional contributions depend on the choice of the 
origin of the coordinate system. This situation is not satisfactory, as a physical observable 
should be origin independent. To rectify this, Neese and coworkers suggested to choose the 
origin differently for each excitation such that these additional contributions are minimized 



58 . Usually, this is equivalent to placing the origin on the atom at which the excitation 
takes place. However, in cases where the dipole intensities are very small compared to the 
quadrupole intensities, this scheme sometimes places the origin far away from the excited 



core orbital 51 , which seems unphysical and affects the resulting intensities significantly. 
Moreover, the scheme will also fail for excitations from core orbitals that are delocalized 
over different atomic centers, a situation which occurs for ligand-edge XAS spectra or for 
metal K-edge spectra in polynuclear transition metal clusters. 

Thus, a theoretical framework for the origin-independent calculation of quadrupole inten- 
sities in X-ray spectroscopy would be desirable. Here, we show that such a formulation 
can be obtained if all contributions to the oscillator strengths that are of the same order 
in the wave vector are included consistently. 

This work is organized as follows. The theory is presented in Section [2] After introducing 



the theoretical framework in Sections 2.1 and 2.2 we revisit the multipole expansion 
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of the transition moments in Section 2J3 Subsequently, in Section 2A this expansion 
is applied for the calculation of the oscillator strengths and we demonstrate that these 
become origin-independent if all terms that are of the same order are included consistently. 
The final equations for the isotropically averaged quadrupole intensities are then derived 



in Section |2.5| This is followed by a description of our implementation of the resulting 
formalism within TD-DFT in Section [3j before we illustrate its usefulness for two test 
cases in Section |4j Finally a summary and concluding remarks are given in Sect. [5] 



2 Theory 

For the theoretical description of spectroscopic processes, quantum chemistry commonly 
employs a semi-classical theory. In this framework, the molecules are described with 
(nonrelativistic) quantum-mechanics, whereas the electromagnetic radiation is treated 
classically (for a discussion, see also Ref. |60]). This theoretical framework is also appro- 
priate for absorption and emission processes in X-ray spectroscopy. Here, we will focus on 
the case of absorption, but the results can be transferred to other types of experiments. 



2.1 Electromagnetic Radiation 

Within the Coulomb gauge (i.e., if one chooses the vector potential such that V ■ A = 0), 
a monochromatic, linearly-polarized electromagnetic wave is defined by the scalar and 



vector potentials 61 63 



4>(r,t) = (1) 
A(r,t) = -A £cos(k-r-ut), (2) 

where the wave vector k points in the direction of propagation and its magnitude is related 
to the wavelength by A = 2n/k, where k = \k\. The angular frequency u is u> = 2nv 



with the frequency u, and frequency and wavelength are related by c = \v = u/k, where 
c is the speed of light. Finally, the polarization vector £ is a real unit vector that is 
perpendicular to the direction of propagation (i.e., £ ■ k — 0). 

From these scalar and vector potentials, one obtains for the electric and magnetic fields, 

E(r,t) = -V<j)(r,t) - \ dA ^ 1 ) = A k£sin(k ■ r - ut) (3) 

c at 

B(r, t) = Vx A(r, t) = A (k x £) sin(fc -r-tut). (4) 

Here and in the following, we are using the Gaussian system of units. The electric and 
magnetic fields are perpendicular to each other and to the direction of propagation and 
are oscillating with angular frequency u and the wavelength A. The amplitudes of the 
electric and magnetic fields are Eq = Bq = Aok. 

The intensity I{ui) of the electromagnetic radiation is defined as the energy flux per area 
through a surface perpendicular to the propagation direction. It can be calculated from 
the Poynting vector [61] , 

S = ^(ExB), (5) 

by taking the absolute value and averaging over one period of the oscillations, 

I( W )= r\S\dt = ^Al = ^Al. (6) 

Jq 07T C 07T 



2.2 Molecules in an Electromagnetic Field 

In the absence of an external electromagnetic field, a molecular system within the Born- 
Oppenheimer approximation is described by the nonrelativistic Hamiltonian 

Ho = J2^T + V (n,...,r N ), (7) 

i=l e 

where the momentum operator is given by p = — iTiV, m e and e are the mass and the 
charge of the electron, respectively, and the potential energy V(ri, . . . , rjy) contains the 



electron-nuclei attraction as well as the electron-electron repulsion. Here and in the 
following, the index i is used to label the electrons. 



An external vector potential can be included in this Hamiltonian via 62 64 
1 



H = y,ir- Pi- -Mn,t) 

^— ' 2m, f 



2m P c 



Pi 



2m P m P c 



A{r h t) Pi 



2m e 6 



J2B(r i ,t)-s i + V(r 1 ,...,r N ) 



;A 2 ( ri ,t) 



2m P c 



J2B(r i ,t)-s i + V(r 1 ,...,r N ) 



where g is the electron ^-factor. In the second line we used that in the Coulomb gauge, 
p i ■ A — A ■ p { . After neglecting the term that is quadratic in A, which is justified for 
weak electromagnetic fields, this can be expressed as 



H = H + U(t), 



where the time-dependent perturbation is given by 



(9) 



U(t) 



eA 
m P c 



2m P c 



- \ c os(fc " r i — wt) {£ ■ pi) — | sin(fc ■ Vi — tut) (k X £) ■ s 



(10) 



Here, we inserted the vector potential and the magnetic field of an electromagnetic wave 
given in Eqs. (|2| and Q. Using sin(x) = — [exp(ix) — exp(— is)], this can be expressed 
in the form 



(j(t) = Uexp(-iut) + U*exp(iut), 



with the time-independent perturbation operator, 



u 



eA 
2m e c 



exp(ifc ■ Vi)(£ ■ p^ + i | exp(ifc • n) (k x £) ■ s, 



(12) 



With this form of the perturbation, we can apply Fermi's golden rule to obtain the tran- 
sition rate (i.e., the rate of change in the probability of finding the molecule in the n-th 
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excited state) [62 63 65 



2yj- „ 2 TfA^ 

r „(u;) = — |(0|[/|n)| 5(u -u 0n ) = ^-^\T 0n \ 2 5{u -u 0n ), (13) 
where we introduced the transition moments 

T 0n = — V <0| exp(ifc • n) (p { ■ S) + i 9 - exp(ifc ■ r<) (k x E) ■ a, In). (14) 

Here, |0) and |n) are the eigenfunctions of the time-independent Hamiltonian ifo with 
Ho\n) = E n \n), and transitions only occur if the frequency of the perturbation matches 
the energy differences between eigenstates of the unperturbed molecule, i.e., for oj = u 0n = 
(E n - E )/h. 

Now, Eq. ^ can be used to eliminate Aq from the equation for the transition rate to 
arrive at, 

4tt 2 

r nM = -^—p |T 0n | 2 /(w) 6(u) - u 0n ). (15) 

The absorption cross section, describing the rate of energy transfer from the electromag- 
netic radiation to the molecule, is defined as 



, T 0n (u)huj An 2 n 2 

cron = / — YT~\ — dtu = ~l=r~ \Ton\ , (16) 



' cE 0n 

where Eq u = E n — Eq. Finally, one usually introduces the dimensionless oscillator 
strengths, 

m eC 2m e 2 

J° n = n 2 2t, a ° n = 2 TP ^ 0ri • 



These are defined as transition rates relative to a harmonic oscillator model 63 65 , which 



fixes the prefactor connecting the absorption cross section and the oscillator strengths. 



2.3 Multipole Expansion 



Calculating the oscillator strengths via the matrix elements of Eq. ( 14 ) would in principle 
be possible, but is cumbersome and in general not feasible. The required integrals are 
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difficult to compute analytically (for a possible approach, see Ref. [66]), and because of 
its dependence on the wave vector k, the operator in To„ is different for each excita- 
tion. Therefore, one usually performs a multipole expansion. The starting point for this 
expansion is a development of the exponential in a Taylor series, 

exp(ifc • n) = 1 + i(k ■ n) - -(k ■ rif + ■■■ (18) 



This is substituted into Eq. (14) and, subsequently, one collects the terms of different 



orders in the wave vector k, i.e., 

T 0n = iff + iff + iff + ■ ■ ■ (19) 

In the following, we will consider terms up to second order in k. Here, \k\ =2tt/X acts as 
the expansion parameter, and we note that for larger wavelengths A, the convergence of the 
Taylor expansion will be faster. For typical molecules and wavelengths in the ultraviolet or 
visible range, the wavelength is large compared to the molecular size, and it is sufficient to 
include only the first (zeroth-order) term in this expansion. This corresponds to assuming 
that the oscillating electric field is constant over the whole molecule. However, for the 
short wavelengths used in hard X-ray spectroscopy this approximation is not adequate 
and higher-order terms need to be included. 

2.3.1 Zeroth order: Electric-dipole moment 

In zeroth order in the wave vector k, we have 

T on = E (°\Pi -e\n) = e- (o|A», (20) 



where we have introduced the electric-dipole moment operator in the velocity representa- 
tion 

A p = — J> (21) 

i 
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By using the relations given in Appendix [XJ the matrix elements of the electric-dipole 
moment operator in the velocity representation can be related to those in the conventional 
length representation as 

(0\fi p \n) = ~^E 0n (0\fi\n), (22) 

where we introduced the electric-dipole moment operator in the length representation 

M = e^rV (23) 



Thus, for the zeroth-order contribution, we arrive at 



rp(0) _ rpfa) _ ; E 0n 
1 0n — 1 0n — 1 * 



(f.(0|£|n>). (24) 



2.3.2 First order: Electric-quadrupole and magnetic-dipole moments 



In the first order in k, we find 

T on = £ E <°K fc • r ')(P« -£)\n) + ^fE <°K fc x £ ) ■ *'!»>■ ( 25 ) 

i i 

The matrix elements in the first term can be split into one term that is symmetric and one 
that is antisymmetric with respect to interchanging the wave vector k and the polarization 
vector £ via 

(0|(fc • ri)(p, • S)\n) =^<0|(fc ■ ri )(pi ■ S) + (k - Pl )(n ■ S)\n) 

+ l(0\(k.r i )(p i .£)-(k.p i )(r i .£)\n). (26) 



From the symmetric first term and using Einstein's convention of implicit summation over 
repeated Greek indices, which we use to label the Cartesian components x, y, and z, we 
obtain 

T on ] = ^^k a £p(0\ri, a Pi,p + Pi, a n,f)\n) = ]- k a 8p (0|Q^|n) (27) 
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where we introduced the electric-quadrupole moment operator in the velocity representa- 
tion 

Qa/3 = — E ('VA/? + Pi,a r i,p) ' ( 28 ) 
e i 

Again, with the help of the relations given in Appendix [A] the matrix elements in the 
velocity representation can be related to those in the conventional length representation, 
and one arrives at 

T on=^k a £ {o\Q a ^\n), (29) 

where 

Qa/3 = e ^ T i,a- r ifi (30) 
i 

is the operator of the electric-quadrupole moment in the length representation. Note that, 



in contrast to most other authors 58,67 , we do not introduce a traceless version of the 
quadrupole tensor here. The traceless definition arises from the expansion of l/\r\ that 
is often introduced in the context of intermolecular interactions, whereas in the case of 



an expansion of the exponential exp(ifc ■ r) considered here the definition of Eq. (30) is 
more natural. Nevertheless, because the wave vector k and the polarization vector £ are 
orthogonal, the diagonal elements of the electric-quadrupole transition moments do not 
enter here and it would, therefore, be possible to alter their trace without consequences. 

For the antisymmetric second term, we can use that (k ■ r^) and {fi i ■ £) commute because 
k and £ are orthogonal and then apply the vector identity 

(a • c)(b • d) - (b ■ c)(a • d) = (a x b)(c x d) (31) 

to obtain, 

i 

= i— V(0|(fc x £)(n x p.) In) = ic(fc x £) ■ (olm'lra). (32) 
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with the (spin-independent) orbital magnetic-dipole moment operator 



m 



2m e c 



( T * X ^ 



(33) 



Thus, this antisymmetric term Tq™ ^ adopts the same form as the last, spin-dependent 



term in Eq. (25) 



(m 3 ) 



On 



izr— x S) ■ (0\gsi\n) = ic (k x £) 



2m P t-? 1 



[0\rh s \n), 



(34) 



with the spin magnetic-dipole operator 



m 



2m P c 



52 qk 



(35) 



Combining the two contributions, we arrive at the magnetic-dipole transition moment, 



Tt ] = i^" J2^ k x £ ) ' <°K r * x Pi) +9Si\n) = ic{k x S) ■ (o|m|n), (36) 

e i 

and the total magnetic-dipole moment operator, 



m 



2m P c 



E 



(r« x pj + 5 ^ 



(37) 



Altogether, the first-order transition moments consist of an electric-quadrupole and a 
magnetic-dipole contribution, i.e., 



T^ 1 ) _ rp{Q) | T (m') ^(m 5 ) _ ^(Q) ^(m) 
J On ~~ 1 On "r J On "r J On — J On "r J On ■ 



(38) 



2.3.3 Second order: Electric-octupole and magnetic-quadrupole moments 



In second order in k, we find, 



n(2) 
On 



(0|(fe ■ rO(fc • rOCPi • f )|n> - ^- ^ (o|(fc • r,)(fc x £) • Si |n>. (39) 



2m e 

In a similar fashion as for the first-order term above, the matrix elements in the first 
term are split into a part that is symmetric and one that is antisymmetric with respect 
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to interchanging the polarization vector £ with one of the wave vectors k, 

(0\{k ■ ri ){k ■ r^p, ■ £)\n) 

= ^{0\(k-r i )(k-r i )(p i -S) 

+ (k ■ n){k ■ ■ £) + (k ■ Pi ){k ■ r^fa ■ S)\n) 



+ l(0\2(k-r i )(k-r i )ip i -S) 



3 



(k ■ r t )(k ■ p^in -£)-{k- p t )(k ■ r,)(r, ■ £)\n). (40) 



For the symmetric first term, we obtain 



= -^E fc ^<0|O^», (41) 

i 

with the operator of the electric-octupole moment in velocity representation 



= {ri,aTi,pPi,i + r^aPi^r^ + Pi, a ri^r i>7 ) . (42) 



Using the relations given in Appendix [A] these matrix elements in the velocity represen- 
tation can be converted to those in the conventional length representation, 

Tq° ] = i^k a k fi £ 7 {o\6 a p 7 \n), (43) 

with the octupole operator in length representation given by 

da?-/ = e ^2 n^n^j. (44) 



Again, note that our definition differs from the one given elsewhere 67 , as we do not 
introduce a traceless form here. In fact, for the Cartesian expansion of the exponential 
exp(ifc • r) it turns out that introducing such a traceless definition here is not possible be- 
cause only the terms depending on the trace of the octupole moments, O aa $ will contribute 
to the isotropically- averaged oscillator strengths later on. In contrast, when describing 
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intermolecular interactions starting from an expansion of l/|f|, these terms are zero and 
do not appear. 

After some algebra (see Appendix [Bj, the antisymmetric part of Eq. (39 ) can be expressed 



as 



(M>) 



On 



6m £ 

e 
6m. 



^2(k x £) ■ (0\(k ■ n) ■ (r< x p { ) + (n x p i )(k ■ Vi)\n) 

i 

^2(k x S) a kp (0|r i)/3 (n x Pi ) a + (n x p^r^n) 



(k x S) a k p (0\M' aP \n), 



with the (spin-independent) orbital magnetic-quadrupole operator 67,68 

M'af, = l ( r iA r i X Pi)* + ( r i X Pi)^,^. 

e i 

Note that this operator is not symmetric with respect to interchanging a and (3. 



(45) 



(46) 



The remaining spin-dependent part of Eq. ( 39 ) is given by 
rriM')__ eg J^ik x £) ■ (0\(k ■ r 2 ) ■ s^n) 



On 



2m P 



i 

with the spin contribution to the magnetic-quadrupole operator 



M 



9 1 n,i3 Si,, 



(47) 



(4* 



2m e c 

Finally, the full magnetic-quadrupole contribution is obtained by adding the orbital and 
spin contributions to obtain 



(M) _ rrW) . ^(M") 



1 0n — 1 0n r 1 < 



On 



-- (k x S) a k p (0\M' a p + M s a p\n). 



(49) 



Altogether, the second-order transition moments consist of an electric-octupole and a 
magnetic-quadrupole contribution, i.e., 



T (2) _ T (0) T (M') T (M°) _ T (0) T (M) 
1 On ~ 1 On ' On ' 1 On ~ 1 On ' 1 On ■ 



(50) 
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2.3.4 Summary of the multipole transition moments 



In summary, the multipole expansion of the transition moments of Eq. ( 14 ) up to second 
order results in five different contributions, 

rp _ rp{fj) , rp(Q) , rp{m) , rp{0) . rp{M) . /r-i\ 

1 0n — I 0n + i 0n + 1 On + 1 On + 1 On l 0i J 

The expressions for these contributions are summarized in Table |T] In zeroth order, one 
encounters the well-known electric-dipole transition moments. Starting from first order, 
an electric and a magnetic contribution appear. In first order, these are the electric- 
quadrupole and magnetic-dipole transition moments, whereas in second order one has the 
electric-octupole and magnetic-quadrupole transition moments. 

[Table 1 about here.] 

If we restrict our considerations to a spin-orbit coupling free framework in the absence of 
static external magnetic fields, the wavefunction can always be chosen as a real function. 
In this case, the electric transition integrals, (0|/t Q |n), (0|Qa,g|n), and (0|O a/ 3 7 |n), as 
defined here are always real. Moreover, for the magnetic transition moments, the spin 
contributions can be neglected because for states with a multiplicity larger than zero 



(i.e., for S > 0) the different Ms-components of the multiplet will be degenerate 69,70 
and the components with +M5 and — Ms provide spin contributions to the magnetic 
transition moments that cancel each other. Therefore, these spin contributions, (0 n) 
and (O ^ n), will not be considered further in the following. The remaining orbital 
contribution to the magnetic transition integrals, (0 \m a n) and (0|vM ai a n), are then 
purely imaginary. As a consequence, we notice that the zeroth-order and second-order 
transition moments are purely imaginary, whereas the first-order transition moments are 
purely real. This holds both for the electric and for the magnetic contributions. 



16 



2.3.5 Origin dependence of the multipole transition moments 



While the electric-dipole transition moments are independent of the choice of the origin 



of the coordinate systems, the higher-order transition moments are origin dependent 67 
We give derivations of the expressions for the change of the electric-quadrupole and electic- 
octupole transition moments and of the magnetic-dipole and magnetic-quadrupole transi- 



tion moments upon shifting the origin from O to O+a in the Supplementary Material 71 
Here, we only present the final results as required for the following discussions. 

For the electric-quadrupole transition moments, one has 

(0\Q af )(O + a)\n) = (0\Q a p(O)\n) - a p (o\fl a \n) - a a {0\fip\n), (52) 

whereas for the electric-octupole transition moments, the corresponding expression be- 
comes, 

(0 + a)\n) = (0\d afh (O)\n) 

- a^(0\Q a p(O)\n) - ap(ti\Q ai {0)\n) - a a (0\Q^(O)\n) 
+ a a ag(o|/i 7 |n) + a a a 7 (o|/tg|n) + a^a 7 (o|/t a n). (53) 

Thus, the electric-quadrupole and octupole moments are only origin independent if all 
lower-order electric transition moments vanish. The above expressions hold both for the 
length and for the velocity representation. 

For the magnetic-dipole transition moment, a shift of the origin of the coordinate systems 
results in, 

(0\m' a (O + a)\n) = (0\m' a {O)\n) - e afil a p ^{o\^\n) 

i E 

= (0\m' a (O)\n) + -e a(Sl a^-^(0\fi^\n), (54) 

where e a/ g 7 is the Levi-Civita tensor. Thus, a term depending on the electric-dipole mo- 
ment in the velocity representation arises. This can be converted into the length represen- 
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tation using the relation from Appendix [AJ The origin dependence of the magnetic-dipole 
transition moments vanishes if the electric-dipole transition moment is zero. 

Finally, for the magnetic-quadrupole transition moment the origin dependence is given 



by 68,72 



(0\M' af3 (O + a)\n) = (0\M' aP (O)\n) 

1 2 

- —e a7S a 1 (0\Qp S (O)\n) + — £ aiS apa 7 (0\pP s \n) 

+ p afS (a • (0|m'(O)|n>) - 2a p (0\m' a (O)\n) . (55) 

Again, the electric-dipole and electric-quadrupole transition moments in the velocity rep- 
resentation can be converted to the length representation to finally arrive at 

(0\M' af} (O + a)\n) = (0\M' aP (O)\n) 

+ g -^7 £ay8a~f{0\Qp5{O)\n) - — — e ajS 0/304 (O^n) 

+ 2 -5 a p (a • (0|ro'(O)|n)) - 2a^\m! a {0)\n). (56) 

Here, we notice that upon shifting the origin, the magnetic-quadrupole transition mo- 
ment generates all lower-order contributions, i.e., terms depending on the electric-dipole 
and electric-quadrupole transition moments as well as on the magnetic-dipole transition 
moments. 



2.4 Oscillator Strengths 

The multipole expansion of the full transition moments T$ n can now be inserted into 



Eq. (17) to obtain an expression for calculating the oscillator strengths, 

_ 2m e ■ (0 ) (1) (2 ) ,2 
JOn — 2T? On ~r 1 On "r J On ~r ' ' ' 

_ 2m e , (At) ( q) (m) (Q) (M) | ,2 

_ 2 On "T" On On On "T" On "T ' ' ' ■ 

e z -&0n 
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Here, different truncations of the expansion can be employed. In the dipole approxima- 
tion, only the zeroth-order term is retained, and for the oscillator strengths one arrives 
at the well-known expression in which the oscillator strengths are proportional to the 
squared absolute value of the electric-dipole transition moments. This approximation is 
commonly employed in electronic spectroscopy in the ultraviolet and visible region of the 
electromagnetic spectrum. 

Here, we are interested in cases where the dipole approximation breaks down, such as 
K-edge XAS spectroscopy of transition metal complexes. In such situations, higher-order 
terms in the multipole expansion have to be included. In the currently used approxima- 
tion, the multipole expansion of the transition moments is truncated after the first order, 



i.e., the oscillator strengths are approximated as 58 



However, it turns out that the resulting expressions depend on the choice of the origin of 



the coordinate system (see Ref. 58 for details). 



[Figure 1 about here.] 



To obtain an origin- independent formulation, we return to Eq. (57) and realize that the 
squared absolute value results in a sum of products of multipole transition moments. 
These products are of different orders in the wave vector k, as is illustrated in Fig. [TJ 
Hence, it seems logical to retain all terms up to a given order in k in the expression for 
the oscillator strengths instead of truncating the multipole expansion of the transition 
moments. By collecting terms that are of the same order, the oscillator strengths can be 
expressed as 

fon = fon + fon + fon H (59) 
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where 



f (o) 

J On 
J On 



2m e 
e 2 E 0n 
2m„ 



(0)1 2 



On 



2Re(T ( ^T «)=0 



(2) 
On 



e 2 E t 
2m e 



On 



e 2 E, 



On L 



T (1) I 2 -I. 9Rp/T w '* 
J 0n I i~ ZXU3 V 0« On J 



(0),* ^(2^ 



(60) 
(61) 
(62) 



where the star denotes complex conjugation. Because Tjj^ is purely imaginary and T ^ is 
real, their product is also purely imaginary and the first-order contribution vanishes. 
In the following, we will retain all terms up to second order and it turns out that the 
resulting approximation for the oscillator strengths is independent of the choice of the 
origin. 



2.4.1 Origin independence of oscillator strengths 



Starting from the definitions of 7q° , T ^% and in Eqs. (20), (25), and (39), respec- 
tively, one can easily see that their origin dependence is given by 



T£\0 + a) = T^\0), 



T^(0 + a) = T^>(0) + i(k.a)TZ'(0 



(o), 



T^\0 + a) = T£(0) + i(fc ■ a) T^(0) - -{k • af T^{0). 



^2^(0), 



(63) 
(64) 
(65) 



Therefore, the zeroth-order contribution to the oscillator strengths , i.e., the expression 
obtained in the dipole approximation, is obviously origin independent. For the second- 
order contribution, we have 

2m P 



f!8(0 + a) 



e 2 E ( 
2m 
^E ( 



(k.a) 2 \T^(0) 



\Tg\0 + a)\ 2 + 2Re(T ( ^*(O) T^(0 + a) 
; T «(O)| 2 + 2Re(i(fc • a)T^(0)T^*(0) 
+ Re(V f *(0) [2T^{0) + 2i(k • a)T«(0) - (k ■ a) 2 T^(0)}) 
f!S(0) (66) 
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and find that this contribution is indeed independent of the choice of the origin. In fact, it 
can be shown that for each order, the higher-order contributions to the oscillator strengths 
are origin independent if all terms that are of the same order in the wave vector k are 
included. This is demonstrated in Appendix [Cj 



2.4.2 Dipole and quadrupole oscillator strengths 

After having established an origin-independent definition of the different approximations 
to the oscillator strengths, we will now turn to deriving explicit expressions. Considering 
only the zeroth-order contribution corresponds to the dipole approximation, in which the 
dipole oscillator strengths are given by 

hn - /i? = f [ oi ] = ^E 0n \E ■ <0|A|n>| 2 = ^ 2 E 0n (S a (0\^\n)) 2 . (67) 

Since the first-order contributions vanish, the next step to go beyond the dipole approx- 
imation is to include all second-order contributions. Thus, the oscillator strengths can 
be approximated as the sum of the dipole (zeroth-order) oscillator strengths and the 
quadrupole (second-order) oscillator strengths, 

hn - /«? + /£> = KJ\ 2 + KJ\ 2 + 2Re(T^). (68) 

We will refer to this approximation as the quadrupole approximation. For the quadrupole 
oscillator strengths, we can insert the individual multipole transition moments, and obtain 
five different terms, 



2m e 



n(2) _ ^""e |t^(Q)| 2 _i_ \rp(m) | 2 , ^-o ( T (Q),*r r (m)\ 9R a {rp(n),*rp(0)\ . 9R a (rp{fi),*rp(M)\ 

JOn ~ 2Z? On | | On \ T^rte^ i Qn ± 0n ) -|- zruj^j 0n i Qn j -t-zne^i 0n ± 0n ) 
e rjQ n l 

— J On ' J On ' J On ~r J On "t" J On ■ \ VJ ) 

First, there are three contributions arising from products of first-order transition moments, 
an electric-quadrupole-electric-quadrupole contribution, 

fSP = ^2 E on <0|g Q/3 |n>) 2 , (70) 
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a magnetic-dipole-magnetic-dipole contribution, 



(m 2 ) 
On 



2m P c 2 



e 2 E t 



On 



(k x £) a (Q\rh a \n) 



2m P c 2 



e 2 E l 



On 



(k x £) a \m(^\m a \n) 



(71) 



and a cross-term, the electric-quadrupole-magnetic-dipole contribution, 

2m e c 



AQm) 
J On 



e 2 h 



k a £(3 (0\Qap\n)j ({k x S) a lm(0| 



m a \n 



(72) 



These three contributions have been considered previously in the calculation of the quadru- 
pole oscillator strengths in Ref. [58). In addition, two additional contributions have to be 
included in order to collect all terms that are of second order and to arrive at an origin- 
independent approximation. These are the electric-dipole-electric-octupole contribution, 

2m P 



J On 



'3h 2 e 2 



E 0n (£ a (0|/t a |n)) (kakpEj (0\O a/ 3 7 \n)), 



(73) 



and the electric-dipole-magnetic-quadrupole contribution, 



f QiM) 

J On 



2m e c 



S a (0\fj, a \n)) ((k x £) a kp lni(p\M a p\n) 



(74) 



Now we choose the wave vector as k = ke x along the x-axis and the polarization vector 
as £ = e y along the y-axis. Consequently, (e x x £) becomes the unit vector e z along the 
z-axis. This is no loss of generality, as the molecule can still have an arbitrary orientation 
in the coordinate system. Using 

Equ 



k 



he 



the different contributions to the oscillator strengths become 



f (M 2 ) 

J On 

AQ 2 ) 

J On 
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2m 



2^0n 



(Qm) 
On 
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(75) 



(76) 
(77) 
(78) 
(79) 



ft 0) = -^ShEl (0\fiy\n)(0\O xxy \n) (80) 



3e 2 h A c 2 ' 
e 2 h 3 c 



tn M) = fS^n (0|^|n)lm(0|^|n>. (81) 



These oscillator strengths refer to an experimental setup in which the incident radia- 
tion has a well-defined polarization and in which the molecules have a fixed orientation 
with respect to the radiation. Using the expressions for the origin dependence of the 



different multipole transition moments given in Section 2.3.5 it can be verified that the 
total second-order oscillator strengths calculated using the above equations are origin- 
independent (see Supplementary Material |71|). 



2.5 Isotropic Averaging 

Often, the molecules are not oriented with respect to the incident radiation in experiments, 
but the measurement is performed in solution where the molecules can freely rotate. Thus, 
to arrive at final expressions for the oscillator strengths in such experiments, we have to 
perform an averaging over all possible orientations of the molecule. 



The expressions for performing this averaging are derived, for instance, in Ref. |67| (see 
in particular chapter 4.2). For the isotropic averages of tensors with two, three, and four 
Cartesian indices, one finds 



(T xx )i SO — 'y^jjigip) iso T a /3 (82) 



a/3 



{T X yz)iso — ^ ^ {iajpk^isaTgpry (83) 



a/37 



(T XX yy^[ so ^ ^ gipj'-fi 5) \so T a /3^$, (84) 
a/37<5 

where the isotropic averages of the Cartesian unit vector i = e x , j = e y , and k = e z are 
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given by 



(iaipj-yjs) 



iso g "af3 i 
1 

iso 7>^a/3 7 ■ 
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(85) 
(86) 
(87) 



For all other tensor components, such as, e.g., (T xy ) iso or (T xxy ) iso , the isotropic averages 



are zero. 



Using these expressions, we obtain for the isotropically averaged electric-dipole-electric- 
dipole contribution to the oscillator strengths, 



(fon ^) 



2m e 
3^ 



2m, 



E 0n 5Z(°l^ a l n ) = ^A jBon ( l^ t l n ) • 

oc n 



(88) 



Similarly, for the electric-quadrupole-electric-quadrupole contribution, we find 



</r> 



ap 7 <) 



m, 



20e^ 4 c 2jBon 



^<o|g Q/3 |n) 2 -^(^<o|g Qa |n) 
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We note that this is identical to the expression in Ref. 58 , where a traceless definition of 



the quadrupole moment is used. For the magnetic-dipole-magnetic-dipole contribution, 
the isotropic average is, 



, ,(m 2 )v 
VOn / 



2m P 



■E 0n ^lm(0|m a |n) 2 = ^^E 0n (lm(o|m|n) 



(90) 



a 

The isotropic average of the electric-quadrupole-magnetic-dipole contribution to the os- 
cillator strengths, 



(fS m) ) 



3e 2 h 3 c 



E 0n ^^a/3 7 (0|Qa/3|^)(0|m 7 |n) = 



(91) 



a/3 7 
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turns out to be zero because (0|Q Q ^p) = (0|(5/3q|^). Finally, for the electric-dipole- 
electric-octupole contribution to the oscillator strengths, we obtain 



(/0n O) )iso = - 45 ^4 c2 ^0n E ( 4 <W<V ~ <WjM ~ $a6$Pj) (0 \p. S \n) (0 1 6 a ^ \ Tl) 



= -£^2 E on E<°I^W(0|Oa^|n>, (92) 

a/3 

where we used the symmetry of the octupole moments with respect to the exchange of 
indices, and for the electric-dipole-magnetic-quadrupole contribution, 

(f^ M) U = ^T c E ln E^(0|^|n>lm<0|^|n>. (93) 

a/37 

Note again that the magnetic-quadrupole transition moments are in general not symmet- 
ric or antisymmetric with respect to the interchange of the two Cartesian indices, i.e., 
(0\M a p\n)^±(0\Mp a \n). 

In summary, there are five contributions to the isotropically averaged oscillator strengths 
up to second order. Also at this stage it can be verified that the resulting total oscillator 
strengths are independent of the choice of the origin, which is shown in the Supplementary 
Material [71] . Note that the individual contributions are still origin dependent. Therefore, 
a separation into electric and magnetic contributions will also depend on the choice of the 
origin. 



3 Computational Methodology and Implementation 

The theory presented here for the origin-independent calculation of quadrupole oscillator 
strengths is applicable in combination with any quantum-chemical method that is capable 
of providing excited states, either via a time-independent formulation or with response 
theory (for a review, see, e.g., Ref. (73])- Here, we select TD-DFT which has become an 
important tool in computational X-ray spectroscopy in the past years [l8|[23|[74|[75] . 
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We have implemented the calculation of the second-order oscillator strengths into the TD- 



DFT module 76 - 78 of the Amsterdam density functional ( Adf) program package 79 80 
Within TD-DFT, the required electric and magnetic transition moments are calculated 
as products of the solution vectors (X + Y) and (X — Y), respectively, with the corre- 
sponding matrix elements in the basis of Kohn-Sham molecular orbitals (for details, see, 



e.g., Refs. 77, 81 821). The required electric-octupole and magnetic-quadrupole integrals 



are provided by Adf's AO Response module 83 



So far, we have only shown that the theory presented here is origin-independent for 
the exact eigenfunctions of Ho. However, an additional difficulty arises in approximate 
calculations. For deriving the equations for the origin-dependence of the magnetic-dipole 



and magnetic-quadrupole transition moments [Eqs. (54) and (56)], we have converted the 
occurring electric-dipole and electric-quadrupole transition moments from the velocity to 
the length representation. However, this is only exact in the case of a complete, infinite 
basis set. Thus, in calculations using a finite basis set, the magnetic-dipole and magnetic- 



quadrupole transition moments do not show the exact origin dependence of Eqs. (54) 



and (56) 



Therefore, we calculate the electric-dipole, electric-quadrupole, and electric-octupole tran- 
sition moments in the velocity representation for the second-order contributions to the 
oscillator strengths (see also Refs [77[[85] for the calculation of transition moments in 
the velocity representation with TD-DFT). It can be easily verified that this results in 
second-order oscillator strengths that are origin-independent also in finite basis-set cal- 
culations. Note that the electric multipole transition moments in the length and in the 
velocity representation are only equal in the basis set limit. However, the calculation of 
higher-order transition moments requires sufficiently large basis sets anyway, so that the 
values in the length and in the velocity representation are usually in very good agreement. 

For the calculation of X-ray absorption spectra in the following, we have employed the 
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scheme of Stener et al. (29) to allow only excitations from the relevant core orbital (see also 



Refs. 30 34 37 for related schemes). For the CI K-edge in TiCU, only excitations from the 
Is orbital of one of the chlorine atoms were considered, while a frozen core was used for the 
other three chlorine atoms in order to obtain a localized core hole |30). For the Fe K-edge 
in vinylferrocene, only excitations from the iron Is orbital were included. All molecular 



structures were optimized using the BP86 exchange-correlation functional 86 , 87 and 
Adf's TZP basis set. The TD-DFT calculations were performed using the BP86 func- 
tional and the TZ2P basis set and employed a fine numerical integration grid (integration 
accuracy 8). All calculations were performed with the scalar- relativistic zeroth-order reg- 



ular approximation (ZORA) [88-91 



4 Results and Discussion 

To illustrate the origin-independent calculation of quadrupole intensities in X-ray absorp- 
tion spectroscopy (XAS) using the theory derived above and to verify our implementation, 
we consider two test cases. As the first example, we use titaniumtetrachloride TiCL- (see 
Fig. |2^i for the molecular structure) and calculate the CI K-edge XAS spectrum. This 



example was considered earlier in Refs. 29,44,581. For such ligand K-edge spectra, the 



prepeak transitions are dipole-allowed, and the second-order contribution to the oscillator 
strength should be small compared to the dipole contribution. 

[Figure 2 about here.] 



For the lowest-energy CI K-edge excitation, the different contributions to the isotrop- 



ically averaged oscillator strengths are calculated using Eqs (88)-(93), and are listed 



in Table [IT] for different choices of the origin. In addition, we included the oscillator 



strengths calculated using the approximation of Ref. 58 , i.e., considering only the electric- 
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dipole-electric-dipole, electric-quadrupole-electric-quadrupole, and the magnetic-dipole- 



magnetic-dipole contributions [cf. Eq. (58)] as well as the full second-order oscillator 
strengths + [cf. Eq. ph]. 



The most natural choice for the origin is the chlorine atom from which the ls-electron is 

(u 2 ) 

excited. In this case, the electric-dipole-electric-dipole contribution to the oscillator 
strength is several orders of magnitude larger than all the second-order contributions, 



and the approximation of Ref. 58 gives results that are identical to the full second- 



order oscillator strengths. In this example, the scheme suggested in Ref. 58 to choose 



the origin such that the sum of the electric-quadrupole-electric-quadrupole /q„ and 
the magnetic-dipole-magnetic-dipole contributions /g n is minimized leads to an almost 
identical choice of the origin. Thus, this scheme is appropriate here. 

The situation changes if the origin is not placed at the chlorine atom. To demonstrate 
this, we moved the origin to the titanium atom. Now, the electric-quadrupole-electric- 
quadrupole ^ and the magnetic-dipole-magnetic-dipole contributions f^™ ^ increase 
significantly and become several times larger than the dipole oscillator strength . 



As a consequence, within the approximation of Ref. 58 the oscillator strength increases 
by more than a factor of two when shifting the origin from the chlorine to the titanium 
atom. However, also the magnitudes of electric-dipole-electric-octupole and the electric- 
dipole-magnetic-quadrupole contributions, and fun M \ increase and since these have 
a negative sign, they exactly cancel the increase of ^ and \ Thus, the full second- 
order oscillator strength remains unchanged. 

In addition, we also shifted the origin away from the molecule by larger amounts. In 
particular, we used shifts of 10 A, 50 A, and 100 A along the negative x-direction. Here, a 
similar observation can be made. The electric-quadrupole-electric-quadrupole, \ and 
the magnetic-dipole-magnetic-dipole, , contributions increase substantially, and for 



a shift of 100 A, the oscillator strength within the approximation of Ref. 58 is four orders 
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of magnitude larger than for the origin at the chlorine atom. On the other hand, when in- 
cluding the electric-dipole-electric-octupole and the electric-dipole-magnetic-quadrupole 
contributions, the full second-order oscillator strengths are unchanged, even though the 
individual contributions differ. 



[Table 2 about here.] 



As a second example, we consider vinylferrocene, which is a ferrocene molecule bearing a 
vinyl substituent at one of the cyclopentadienyl rings (see Fig. for the molecular struc- 
ture). Here, we consider the Fe K-edge XAS spectrum and specifically the lowest-energy 
(prepeak) excitation, which is a Is — > 3d transition. In unsubstituted ferrocene, this 
prepeak excitation is dipole-forbidden for symmetry reasons, and its oscillator strength 
is solely due to the second-order contributions. In this case, the electric-quadrupole- 
electric-quadrupole and the magnetic-dipole-magnetic-dipole contributions become origin 



independent (see Section 2.3.5), whereas the remaining second-order contributions 



I On 



and fon M ^ vanish. However, in vinylferrocene this symmetry is lost and the lowest-energy 



transition gains a small dipole oscillator strength (for a detailed discussion, see Ref. |5T]). 
The oscillator strengths and their contributions calculated for the lowest-energy Fe K-edge 



excitation using different choices of the origin are shown in Table III First, the most nat- 
ural choice for the origin is the iron atom. In this case, the electric-dipole-electric-dipole 
and the electric-quadrupole-electric-quadrupole contributions to the oscillator strength 
are comparable in size. The remaining contributions are orders of magnitude smaller. 



Therefore, the oscillator strength calculated with the approximation of Ref. 58 is iden- 
tical to the full second-order oscillator strength. 

To investigate the dependence on the origin, we shifted the origin far away from the 
molecule using a shift of 100 A in the negative x-direction, a shift of 100 A in the neg- 
ative z-direction, and a shift of 50 A in both the negative x-direction and the nega- 
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tive z-direction. In all three cases, the electric-quadrupole-electric-quadrupole and the 
magnetic-dipole-magnetic-dipole contributions, ■* and /q™ ^ , increase by several orders 
of magnitude compared to the calculation in which the origin is placed at the iron atom. 
As a result, the oscillator strengths calculated with the approximation of Ref. |58| also 
increase by up to five orders of magnitude. However, at the same time the two remaining 
second-order contributions, i.e., the electric-dipole-electric-octupole contribution /g^°' > 
and the electric-dipole-magnetic-quadrupole contribution /on^, assume large negative 
values and exactly cancel the increase of /g„ and / 0n such that the total second-order 
oscillator strength remains origin independent. 

Finally, we used the scheme suggested in Ref. |58| for fixing the origin of the coordinate 
system, i.e., we chose the origin such that the sum of ^ and /q™ ^ is minimized. In 
the situation considered here, where the electric-dipole-electric-dipole and the electric- 
quadrupole-electric-quadrupole contributions to the oscillator strengths are of similar 
size, this scheme moves the origin away from the iron atom. The resulting shift is given in 



the caption of the last column of Table III As a consequence, the oscillator strength within 



the approximation of Ref. 58 decreases by ca. 30 %. Again, this decrease is compensated 



if the remaining second-order contributions are included. Thus, the scheme of Ref. 58 
can lead to a spurious decrease of the oscillator strength in some cases. Previously, we 
found that this problem is even more severe in cases where the electric-dipole-electric- 
dipole contribution to the oscillator strength is significantly smaller than the quadrupole 



oscillator strength 51 . However, if all second-order terms are included consistently the 
quadrupole oscillator strengths become origin-independent and no special placement of 
the origin is necessary. 



[Table 3 about here.] 
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5 Conclusions 



We have derived origin-independent expressions for calculating XAS intensities beyond 
the dipole approximation. In particular, we have shown that for a consistent formulation, 
it is necessary to retain all contributions to the oscillator strengths that are of the same 



order in the wave vector. This differs from the previous approach 58 , in which the 
multipole expansion was truncated for the transition moments. Here, two additional 
contributions to the second-order (quadrupole) oscillator strengths arise, which are cross- 
terms depending on products of electric-dipole and electric-octupole transition moments 
and of electric-dipole and magnetic-quadrupole transition moments, respectively. 

Thus, the origin dependence of the sum of electric-quadrupole-electric-quadrupole and 
magnetic-dipole-magnetic-dipole contributions pointed out earlier |58| is not a funda- 
mental limitation of the use of the multipole expansion. In fact, we could show that 
to arbitrary order in the wave vector, origin-independent expressions for the oscillator 
strengths are obtained if all terms of the same order are included consistently. Conse- 
quently, within the multipole expansion it should always be possible to derive origin- 
independent expressions for physical observables. 

An origin- independent formalism for calculating quadrupole intensities is particularly im- 
portant for studying ligand and metal K-edge XAS spectra of transition metal complexes. 
To this end, we have implemented our theory for calculating XAS spectra with TD-DFT, 
and applied it to two simple test cases. Here, we want to stress that our results do not 



invalidate any previous results obtained with the formalism of Ref. 58 . On the contrary, 
our test calculations showed that the two additional contributions are negligible as long 
as the origin of the coordinate system is placed at the atom where the core excitation oc- 
curs. However, with our origin-independent theory, it is no longer necessary to make sure 
that the origin is chosen appropriately. This is particularly important for cases where the 
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quadrupole intensity is larger than or comparable to the dipole contribution, where the 
scheme proposed in Ref. [58] might place the origin far away from the relevant core orbital. 
Moreover, it makes it possible to treat excitations from core orbitals that are delocalized 
over several atoms (e.g., for calculating ligand K-edge spectra or metal K-edge spectra in 
polynuclear transition metal complexes) without the need to perform a transformation to 
localized core orbitals. 

Of course, the theory presented here is not limited to TD-DFT, but can be employed 
for the calculation of quadrupole intensities in combination with any quantum-chemical 
method capable of providing the required transition moments. Moreover, it is not re- 
stricted to XAS spectroscopy, but is also applicable for calculating XES intensities, for 



instance using the approach of Ref. 59 . Finally, we note that it becomes necessary to go 
beyond the dipole approximation, not only for short wavelengths, such as those employed 
in hard X-ray spectroscopy, but also for extended molecular systems. For describing 
the optical response of an extended nanostructure in the visible spectrum, it becomes 
necessary to go beyond the dipole approximation as well. Thus, the origin-independent 
formalism derived here will also be essential for predicting optical properties of nanos- 



tructured materials, such as, for instance, metamaterials 92 . 
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Appendix 



A Length and velocity representation 

To show how the electric-multipole moments in the velocity representation can be con- 
verted to those in the conventional length representation, we use the following commuta- 
tors of (products of) the Cartesian components of the position operator with the molecular 
Hamiltonian given in Eq. ([7]), 

[r i;a ,H ] = — p ia , (A-l) 
m 

ih 

[ri^r^p, H ] = —(Pi, a r it p + r ija p i<l3 ) , (A-2) 
[ri,an,i3ri,y, H ] = —{f>i, a r ifi r in + r^ a p ifi r ia + r^r^p^) . (A-3) 

Next, we employ that the matrix elements of the commutator of an operator A and H 
are given by 

(0|[i,iio]|n) = (0\AH - H A\n) = E n (o\A\n) - E (o\A\n) = Eo n (o\A\n) (A-4) 

Here, it is important to point out that this relation is only valid for the exact eigenfunctions 
of Ho and that it only holds approximately for approximate wavef unctions. 

Now, we can use these results to obtain 

(0|pi, a |n) = — (0|[r i>a ,i?o]p) = -iEon-^(0\r i)a \n) (A-5) 
and get for the electric-dipole transition moments 

(Q\fi p a \n) = —(0\p i>a \n) = -i-^-e'Y] (0\r itOl \ii) = -i-|^(o|/t a |n). (A-6) 

i i 

Similarly, we find for the electric-quadrupole transition moments 

(u\QU n ) = - [ ^(°\Q^\ n ) (A-7) 
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and for the electric-octupole transition moments 

(0|O^» = -i^<0|O^». (A-8) 



B Antisymmetric second-order term 



First, the matrix elements in the antisymmetric term in Eq. (40), can be split as 



2{k ■ n)(k ■ r t )(p t -£)-{k- n){k ■ p t )(r t ■ £) - (k ■ Pi ){k ■ r t )(r t ■ £) 

= (k ■ r t ){k ■ ri )( Pi -£)-{k- r t ){k ■ p t )(n ■ £) 
+ (k ■ n)(k ■ rO(p< -£)-{k- Pi){k ■ r^n ■ £). (B-i) 

For both terms, we can employ that both (k ■ r*j) and (p i ■ £) commute because k and £ 



are orthogonal, and subsequently use the vector identity of Eq. (31). In the same fashion 



as for the antisymmetric first-order contribution, we then obtain for the first term, 



(k-r i )[(k-r i )ip t -S)-(k-P i )(r i -e) 

= (k ■ n){k x £) ■ (n x Pi ) = (kx£)-(k- n) (r< x Pi ), (B-2) 

and for the second term, we get, 

(k'r i )(k-r i )(p i -e)-(k-p i )(k'r i )(r i 'e) 

= (k ■ r 4 )(fc • £){k ■ n) - (k • p,)(r, - £){k ■ r,) 

= [{k • n)^ • B) - {k • pt)(n • B)]{k • n) 

= (kx£)-(r i xp i )(k-r i ). (B-3) 
Altogether, we arrive at 

2(k • r 4 )(fc ■ r,)(fr -£)-(k- r,)(* • p 4 )(r 4 • £) - (k ■ p t )(k ■ r t )(r t ■ £) 

= (kx£)- [(k ■ n) • (r< x Pi ) + (n x Pi )(k ■ n)] . (B-4) 
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C Origin independence in arbitrary order 



In Section 2.4.1 , we showed explicitly that the second-order oscillator strengths are inde- 
pendent of the choice of the origin. Here, we prove that this still holds for an arbitrary 



order. From the definition of the full transition moments [Eq. (14)], we find for its change 
upon shifting the origin from O to O + a, 



T(0 + a)= exp(ifc • a)T(O) = ( £ • a)») ( £ r<»>(0; 

n=0 ' n=0 

and can identify the terms that are of order m in the wave vector, 

m .„ 

T (m) (0 + a) = ( fc • a )" T (m - n) (0). 



(c-i) 



n=0 



(C-2) 



For the oscillator strengths, the terms that are of order m in the wave vector are, 
2m P 



f (m \0) 



e 2 E 0n 



^T {n \0)[T {m - n \0 



(C-3) 



n=0 



and when the origin of the coordinate system is shifted, this becomes 
f m) (O + a) = T(n) (° + a ) [T {m - n) (O + a)} * 



On 



2m P 



e 2 E { 



o» 



n=0 



n=0 



p=0 



'-i V? 



J2 -~r( k ■ a ) q [r {m ~ n - q \o) 



g=0 



2m P 



m n m—n 



e 2 E t 



Oil 



( -1 ) 9 T7 ( k - a ) P+q T (n " p) (0) [T (m - n " 9) (0) 



n=0 p=0 g=0 



(C-4) 



Now we eliminate p and g by introducing the new indices r = p + q and s = n — p to 
arrive at 



f {m \0 + a) 



2m P 



m n m—s 

EEE 



' j \r— n+s jr 



un n=0 s=0 r=n-s v ' v ' 



(fc-a) r T (s) (0) [T {m - r " s) (0) 



(C-5) 
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The three sums can be rewritten and put in a different order, which leads to 

„ m m—r s+r j 



e 2 E 0n „ „ J ^ (n - s)\(r - n + s)\ 

■r=0 s=0 n=s v ' v ' 



o m—r r , -. r _j . r 

un r=l s=0 t=0 v ' 



(C-6) 



That the summation here is equivalent to the one in Eq. (C-5) can be seen easily by 
considering the six inequalities corresponding to the sums in the two cases and showing 
that these are equivalent. In the second line above, we have taken the term m = out 
for the first sum and introduced the new index t = n — s. 

Finally, we can use the binomial theorem to realize that, 

that is, the last term in the above equation is zero and, thus, we have shown that in any 
order k, the oscillator strengths are origin-independent. 
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Figure 1: Schematic illustration of the different terms arising from the squared absolute 



value in Eq. (57). The entries in the table indicate the order of the different terms in the 



wave vector k. We retain all terms up to second order, as indicated by the red line. 
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Figure 2: Molecular structures of the model systems considered for the calculation of 
X-ray absorption spectra, (a) Titaniumtetrachloride (TiCL*) and (b) Vinylferrocene. The 
orientation of the molecules within the coordinate system is also indicated. 
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Table I: Overview of the different contributions appearing in the multipole expansion of 
the transition moments of Eq. ( 14 ) up to second order. 
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Table II: X-ray absorption oscillator strength for TiCl 4 , calculated for the lowest-energy 
transition at the CI K-edge (excitation energy 2755.6 eV). The total isotropically averaged 
oscillator strength and its different contributions are given for different positions of the 
origin of the coordinate system. 
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Table III: X-ray absorption oscillator strength for Vinylferrocene, calculated for the lowest- 
energy transition at the Fe K-edge (excitation energy 7051.3 eV). The total isotropically 
averaged oscillator strength and its different contributions are given for different positions 
of the origin of the coordinate system. 
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